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A  Comprehensive  Analysis  Algorithm 
for  Ablation,  Temperature  Fields,  and 
Thermal  Stresses  of  Carbon-based  Noset  ip  Materials 

Huang  Zhenzhong 


ABSTRACT 


•  /  /  *4~  /  f 


This  paper  narrates  the  compos i t e  calculative  principles  f  or  -abl  a t  i  on  , 
conduction,  and  thermal  stresses  of  flight  vehicles  in  a  re-entry 
environment,  as  well  as  a  few  techniques  beyond  these  principles.  In  regards 
to  the  chemical  reactions  of  carbon-based  material  ablative  surfaces,  a 
selective  calculation  result  is  used  to  simplify  the  thermochemical  ablation 
calculations.  A  spaced  and  non -spaced  three-level  explicit  difference  scheme 
was  chosen  for  use  in  the  calculation  of  temperature  fields.  Thermal 
stresses  were  calculated  using  the  Jr  i  n  i  te  El  emen  t  Method  within  the 
Hypothesis  of  Linear  Elasticity,  based  on  the  i ncomp 1 e te  Pt i nc i p 1 es  of 
Variational  Potential  Energy,  wherein  it  is  possible  to  calculate  the 
influences  of  many  different  types  of  displacement  boundary  conditions. 

During  the  compilation  of  the:  compr  ehens  i  ve  analysis  procedures^  here  i  n  ,  an 
automated  mesh  dividing  technique  was  found  which  has  application  in 
selective  calculations.  ' 


'i.~  document  received  November  12,  1^32. 


I  .  Foreword 


The  use  of  carbon-based  materials  to  serve  as  heat  shields  in  fl i oh t 
vehicles  has  long  been  a  primary  objective  of  research  [1-3].  One  of  the 
leading  reasons  for  this  is  to  make  thermal  ablation  an  insignificant  factor 
of  external  change  that  takes  place  during  the  time  of  re-entry.  However, 
the  abrupt  thermal  environment  and  subsequent  onslaught  of  thermal  stressors 
is  a  very  serious  danger.  Resultingly,  a  comprehensive  analysis  algorithm 
that  simultaneously  calculates  thermal  environment,  ablation,  conduction,  and 
thermal  stress  has  been  created,  and,  in  regards  to  the  problems  of  combating 
heat  with  car bon -based  noset  ip  materials,  has  become  a  vital  necessity.  This 
paper  recounts  the  calculation  methodology,  and  although  it  is  equally 
applicable  for  use  with  QIAG-type  nosetips,  the  figures  and  examples  provided 
are  for  the  SAI -type  nosetips  CTN:  transliterations  for  QIAO  and  SAI  may  also 
be  pronounced  as  KE  and  SE  respectively;  literal  translations  for  the 
roman izat ion  are  jacket/case/coated  for  QIAO  and  shell  for  SAI]. 

The  thermal  stress  state  of  heat  resistant  nosetips  not  only  impinges 
upon  the  nosetip's  construction  and  form,  but  also  impinges  upon  the  thermal 
environment,  ablation,  and  changes  that  take  place  within  an  instantaneous 
heat  field.  These  factors  each  have  a  mutual  influence.  Herein,  we  endeavor 
to  fully  link  these  factors  together  and  realize  a  possible  simulation  of 


actual  conditions. 
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II.  Ablation  and  Ablative  Contour 


During  the  course  of  re-entry,  the  carbon-based  protective  heat  shield's 
thermochemical  ablation  primarily  occurs  at  extreme  surface  temperatures  that 
surpass  3000  degrees  K.  During  this  time,  the  carbon  surface's 
thermochemical  reaction  can  be  expected  to  depart  from  normal  theoretical 
conditions  14]  C  TN :  sic].  The  surface  temperature  of  the  heat  resistant 
shield  rises  extremely  fast  during  the  time  of  re-entry,  although  during 
initial  re-entry  it  is  influenced  by  abnormal  thermochemical  reactions. 
However,  for  quantitative  comparisons  of  total  ablation,  this  particular- 
stage  of  ablation  is  very  limited.  Nonetheless,  during  the  time  of  actual 
re-entry,  all  calculations  of  ablation  are  in  accordance  with  normal  chemical 
hypotheses  and  could  not  lead  to  errors  that  would  influence  outcome  of 
calculation  CTN:  sic]. 


The  meticulous  analysis  conducted  for  thermochemical  ablation  of 
carbon-based  noset i ps  uses  thermodynamics  data  taken  from  bibliographical 
reference  C5]  CTN:  the  document  cites  reference  5  when,  in  fact,  it  is 
referring  to  reference  6!].  In  regards  to  the  thermodynamics  numeric  data 
provided  for  five  types  of  vapors  and  10  types  of  compounds:  C  (gaseous) ,  CA , 
CJt  Cr  Cj.,  CO,  C0Z  ,  C^O,  CN,  CN^C  C-N-N)  ,  CNX(N-C-N>,  C^N,  C^N^  ,  and  Cy 

N  ,  a  dressing  by  screening  calculation  was  conducted  for  actual  ablative 


surface  edges  of  the  projectile.  The  results  can  be  seen  in  Figure  1.  Here, 


15)  expresses  successive 


aforementioned  compound  reactions.  Carbonization 
only  a  fraction  of  the  total  thermodynamic  ablati 


formation  of  the 
in  ablative  change  occupies 
ve  rate.  Calculations 


5 


clearly  indicate  that  in  an  ablative  carbon  surface,  the  formation  of  CO. 


C  N ,  and  ON  in  a  reaction  (cor re 


spond i ng  to  .if.  an  d  f?  .  c 


curve  i n 


figure  1  '  ,  in  regards  to  the  carbon's  ablation  rate  or  the  contributions  o 
the  thermal  effect  when  compared  to  other  reactions,  is  much  larger.  This 
fact  alone  could  greatly  simplify  the  calculation  of  ablation.  It  is 
sufficient  for  us  to  calculate  these  major  reactions  and  corresponding 
chemical  elements,  so 


30 

C*K' 

N* 

2CtirV 

As  the  thermodynamic  data  clearly  shows,  in  common  ablative  surface 
temperature  ranges,  the  normal  constants  and  thermal  effects  of  these 
reactions  can  be  completely  expressed  in  Figure  1  with  use  of  similar 
equations  and  approximation  averages  used  as  their  expressions.  This  make 
calculation  of  ablation  much  more  convenient.  The  normal  constant  in  a 
car  bon -oxygen  reaction  can  only  use  an  oxygen  residual  that  is  equal  to  re 
as  a  substitute.  This  is  because  from  the  normal  constant  we  know  that 
oxygen  in  an  ablative  reaction  has  an  actual  residual  value  that  is  so  sma 


SCO 

*CN 


a  a 


we  can  consider  it  as  negligible  from  the  start. 


** - 5 - 5 - T  3  5  J  meters 

Figure  1.  Various  carbon ic  ablative  rates  for  reactions  at  ablative  wall 
represented  as  a  fraction  of  projectile  change. 

The  ablati  ~  a  '  1  •'  s  quantitative  thermochemical  ablative  rate  can  be 
expressed  as 


<pqa&/h. 


where  i  n 


the  dimensionless  ablative 


rate 


3 


derived  as  a  normal  constant  b 


means  of  a  co-routine  to  -figure  the  dispersion  equation  tor  the  attached 


surface  layer. 


The  thermal  current  entering  the  ablative  wall  is 


The  second  and  third  orders  of  the  equation  on  the  left  represent  the 
atmosphere  opposing  the  ablative  wall  and  the  outer  boundary  radiant  heat  o-f 
the  ablative  wall.  Furthermore, 


We  need  only  provide  calculations  -for  denudation  and  particle  immersion, 
then  we  can  per -form  a  comp  rehens  i  ve  calculation  sequence  -for  conduction  and 
contour  changes  [  TN :  sic]  as  influences  of  thermal  stress. 

The  noset  ip,  during  the  course  of  the  ablative  process,  undergoes  an 
instant  exterior  contour  change,  and  this  can  be  represented  as  an  equation 

dr,  _  _  A  rfim 

dt  hi  [  rmr I  p  4 ) 

Here,  is  the  quantative  rate  for  ablation,  denudation,  and  immersion  C  TN : 

particle  immersion].  This  is  a  nonlinear,  hyperbolic  differential  equation. 
The  methodology  used  in  deriving  this  equation  is  first  to  begin  with  an 
ao:  mat i on  linearization  equation,  then  use  the  Implicit  Difference 

Scheme.  This  serves  as  a  type  of  man-made  patch  that  allows  us  to  create  a 
difference  equation.  Finally,  we  use  a  method  of  substitution  to  derive  our 
solution.  B i bl i ogr aph ical  references  E5,  7]  iterate  the  p: : .  ss  in  precise 
detail.  In  Figure  4  at  the  back  of  this  paper,  a  simultaneous  calculation  of 
temperature  fields  is  ou t 1 ined.  These  methods  [TN:  refers  to  above  methods] 
are  used  to  calculate  the  course  of  ablative  heat  protection  for  the  noset  ip 
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III. 


Instantaneous  Temperature  Fields 


The  majority  of  carbon-based  materials  assume  an  an  i  stroph  i  cal  form,  and 
as  a  result  o-f  this  we  must  conduct  calculations  by  means  of  a  thermal 
conduction  equation  for  the  anistrophic  materials.  In  calculating  the 
temperature  fields  for  the  SAI -fabr i rated  nosetip  model,  we  must  break  down 
the  calculations  for  the  3AI-model  nosetip  proper  and  the  secondary 
components.  For  this,  we  will  choose  to  use  an  arbitrary  boundary  surface. 
For  example,  a  spherical  surface  which  meets  a  man-made  boundary  as  depicted 
in  Figure  2.  In  regards  to  the  main  component,  a  spherical  coordinate  system 
was  selected  with  the  corresponding  coordinate  origins.  Furthermore,  a 
mobile  coordinate  system  was  used  in  conjunction  with  the  ablation 
displacement  boundaries  so  the  below  illustrated  coordinate  transformations 
could  be  carried  out: 


With  the  ablation  displacement  boundary  transformation  equations,  we  now  have 
a  stable  boundary.  Furthermore,  thermal  conductivity  for  the  three  mutually 
perpendicular  primary  angles  is  represented  in  the  equation  as 


10 


which  becomes 


iCpCH-llA  d(pCT ) 
A  dt  d£ 


1  3 

sin  $  dO 


(rAsn\9q,)  + 


A  ^  (r***)] 


Here,  i  s  the  rate  at  which  the  width  changes.  This  is  from  the 

external  change  rate  3rs  /  (?*  ■from  equation  (4).  Now  the  ablation 
boundary's  thermal  conductivity  conditions  are 


£-i. 


qt=M~Q 


<  8 


and 


Q-«?.-+p([r“cdryu-ci+c;)-i^.  ■ 

r#  /  3/ 


Now  we  need  only  check  the  zero  angle  of  attack  in  the  body  of 
revolution.  Furthermore,  at  the  very  minimum,  we  must  also  check  the  wind 
axis  and  the  anistrophic  similarities  of  the  level  surface.  With  the 
mater i al s  of  k*  =  k^  ,  we  can  cause  a  differential  abatement  involving^  and 
simplify  the  problem  of  symmetrical  wind  axes.  In  this  type  of  problem  with 
wind  axis  symmetry,  after  the  change  has  occurred  in  equation  (7),  the 
derivatives  of  the  temper ature-to-space  coordinates  of  the  first  and  second 
order,  and  the  thermal  flow  opposing  these  atmospheric  coordinates,  all  have 
internal  CTNs  reference]  points  that  have  a  second  order  precision  in  a 
centered  difference  equation.  We  can  have  a  second  order  precision  for  both 


pre  and  post  differences. 


In  regards  to  the  derivative  of  time,  the 


temperature  field  is  obtained  throught  use  of  the  DuFor t-Franke 1  difference 
scheme.  Afterwords,  from  the  equation  we  get 

ni-^cnv+rm+ocAi1) 

C  <  1  0  > 

Here,  superscript  n  represents  the  time  gap.  Omitting  the  minute 
quantities  of  higher  order,  we  can  correctly  obtain  a  direct  and  simple 
calculation  of  stability  using  the  three-level  explicit  difference  equation. 
To  follow  along  with  the  detailed  process  and  reasoning  of  this  equation, 
refer  to  b  i  bl  i ograph i cal  reference  C  5 1 .  Therein,  it  provides  an  explanation 
on  the  precision  and  analyses  for  this,  the  algorithm  which  provides  us  with 
an  excellent  degree  of  precision  and  influences  the  time  gap  very  little. 


In  regards  to  the  SAI -model  component,  it  is  nonexistent  in  the 
displacement  boundary.  This  is  because  a  fixed  columnar  coordinate  system 
was  chosen  for  use.  We  also  have  a  zero  angle  of  attack  in  the  body  of 
revolution,  so  for  k<^  =  k^  ,  the  thermal  con due*  i ty  equation  can  be  written 
in  ax i symme tr i cal  form  as 


KpCT)  _ 

dt 


r. 


-U  rx  *) 

3r7\  '  dr,  ) 


for  the  main  body  of  the  SAI  component.  For  the  man-made  boundaries  of  the 
SAI  component,  we  must  use  a  method  of  calculation  that  coordinates  thermal 
flow  calculations  with  the  sequence  of  conditions  that  arise. 

In  carrying  out  our  data  analyses,  we  should  point  out  a  few  facts  for 
the  use  of  the  difference  equation  (11):  speaking  in  terms  of  the  $ 


coordinates  within  the  spherical  coordinate  system,  the  main  component  is 


partitioned  into  mesh  lines  based  on  averages 


However,  this  mesh  falls 


within  the  SAI  component's  columnar  coordinate  system  and  is  represented  as 
orientation  and  its  values  are  not  based  on  averages.  Furthermore,  for 
the  secondary  boundary  readings  of  thermal  flow,  this  method  sequenced  the 
conditions  in  the  proper  order.  The  mesh  points  partitioned  within  the  SAI 
component  fall  within  this  boundary  and  are,  of  course',  calculated  with  ths 
main  component  mesh  points.  This  forced  us  to  conform  these  mesh  points  tc 
the  non -spaced  mesh  points  of  the  SAI  component.  At  the  same  time,  the 
boundary  itself  was  also  considered  irregular.  In  order  to  coordinate  this 
with  the  main  component,  for  these  types  of  non-spaced  mesh  points,  a 
three-level  explicit  difference  scheme  must  also  be  used.  It  is  done  by 
performing  the  Taylor  spread  on  the  origin  node's  temperature  and  thermal 
flow.  By  doing  it  in  this  manner,  we  get  a  derivative  for  temperature  and 
thermal  flow  of  non-spaced  gaps  for  a  difference  approximation.  For  exampl 
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Finally,  we  use  a  similar  me  thodol  ogy  for  the  non -spaced  mesh  points  of  the 
main  component  so  we  can  get  a  non -spaced  mesh  point  to  be  used  in  the 
three-level  explicit  difference  scheme. 

-  n,.  +  k,  f(T  (J5) 

or,  02 

For  the  boundary  conditions  found  at  the  irregular  boundary,  a  very 
complicated  difference  scheme  was  constructed.  This  was  due  to  the  linear- 
differences  between  mesh  points  and  the  boundary.  Figure  3  provides  an 
illustration  of  the  boundary  mesh  points.  We  need  a  non -spaced  difference 
expression  for  the  temperature  at  node  <J,1).  Using  a  partial  differential 
of  the  first  order  with  second  order  precision,  we  must,  at  the  same  time, 
also  have  provisions  for  the  intersection  of  the  second  order  partial 
differential  in  the  non-spaced  difference  equation.  This  establishes  the 
difference  equation  for  boundary  conditions  found  in  equation  <15>.  Due  to 
the  use  of  a  non-spaced  difference  approx imat i on  in  the  partial  differential 
for  the  boundary,  aside  from  the  difference  approximation  differential 
itself,  we  must  not  rely  too  heavily  on  man-made  approximations  or  else  our 
results  will  not  as  accurate  CTN:  sic]. 


Figure  3.  Mesh  points  <J,1)  within  the  Mesh  Point  System, 
a :  n/'a 

It  is  in  this  manner  that  we  realize  a  comprehensive  analysis 
algor  it*-.  The  calculations  shown  are  a  model  for  spaced  and  non-spaced 
difference  schemes.  The  time  interval,  stability,  and  precision  are  of 
little  influence  under  these  ordered  conditions. 


Figure  4  illustrates  the  calculation  results  tor  re-entry  ablation  and 
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Figure  4 . 
conduc  t i on 


on  in  a  SAI-mode!  graphite  noset  ip. 


Cal  cu-1  at  i  on  results  -for  re-entry  ablation  and  thermal 
at  varying  altitudes  -for  the  SAI -model  graphite  noset  ip. 
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Car  bon -type  materials  have  a  very  small  range  of  plasticity.  Research 
in  this  area  has  been  limited  to  problems  dealing  with'  strain  and  linear- 
elasticity.  In  the  -fundamental  calculations  for  known  temperature  -fields, 
stress  -fields,  and  strain  -fields,  this  has  been  based  on  the  incomplete 
Generalized  Potential  Energy  Theory  [33  and  it  fits  into  boundary  load  ana 
boundary  displacement  conditions  only  in  terms  of  the  Principles  of  Linear 
Elasticity  as  related  to  minimal  displacement 

n  =  j  j  j [A{en)~ FiUiJdv—  j  j ^Pu,ds 

s; 

v.  1  & 

_ j  |(“i  —  Bf)a,,mds 

Here,  the  Finite  Element  Method  was  used  to  carry  out  the  calculations. 
Within  this  equation,  the  subscripts  i ,j  express  the  weight  and  summation 
sign.  In  problems  of  symmetrical  axes,  we  can  simplify  values  within  the 
columnar  coordinate  system  so  our  problem  becomes  one  of  a  meridian  surface. 
In  the  selection  of  a  triangular  form  to  become  a  finite  element,  aside  from 
d  i  sp  1  acemen  t  If"  at  the  axis  orientation,  the  nodes  are  designated  as 
constants.  This  causes  the  ablative  surface  shield  to  serve  as  the  odd 
integral  in  the  symmetrical  axis.  According  to  the  proposals  cited  in 
bibliographical  reference  [9],  we  can  directly  obtain  values  for 
circumferential  stress  by 


Within  the  singular  element,  we  use  the  fabricated  linearization  to  derive  a 
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Here,  the  superscript  g.  expresses  a  value  -for  the  singular  element, 
the  ordered  nodal  constant  previously  ordered  as 

{<3}*=|,»1  Oi  -js,  v,  wm  vmyT 


We  use  the  extracted  value  -function  -from  equation  (  1 8 >  so  we  may  have  an 
integral  tor  equation  (16).  Next,  a  similar  extraction  is  done  -for  the  other- 
elements  -forming  a  composite,  overall  -functional.  With  an  obtained  numeric 
value  equal  to  double  the  nodal  totals,  the  unknown  quantity  is  spec i tied  as 
nodal  constants  IS  and  u)  which  also  gives  a  sub-equation  -for  the  symmetrical 
coefficient  gap  of  linear  elasticity.  Based  on  this,  we  can  obtain  values 
for  strain  and  stress  from  the  strain-displacement  and  stress-strain 
derivatives. 

In  equation  (16),  because  the  boundary  conditions  are  fully  encompassed 
within  the  functional,  this  in  no  way  causes  illogical  c al cu 1  a t i on  of  those 
boundary  conditions.  Moreover,  we  are  still  in  the  process  of  researching 
the  effects  produced  for  many  different  types  and  kinds  of  boundary 
conditions.  In  this  aspect,  we  can  only  rationalize  the  sustained 
atmospheric  pressure  effects  for  the  externa)  boundary.  However,  the  most 
important  of  all  considerations,  that  of  free  boundary  conditions,  cannot  be 
rationalized.  This  is  a  problem  for  which  we  must  find  a  solution.  In 
coming  up  with  a  proposal  that  meets  all  boundary  conditions,  it  could 
possibly  be  broken  into  three  conditions  or  categories  for  such  an  algorithm. 

The  first  would  be  for  free  boundary  conditions;  the  second  would  be  for  ere 
and  post  boundary  displacement,  for  example,  dealing  with  pre-defined 
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f  or 


aisoUcement  limit  values  of  zero;  the  third  could  be  -for  partial  limits, 
example,  slippage  angle  loss  allowances,  exclusion  of  orientation 
displacement,  etc.  With  these  foundations,  we  can  also,  of  course,  very 
easily  rationalize  the  co-existence  or  unity  of  the  three  afore mentioned 
boundary  conditions. 
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U .  Selective  Calculation  Methodology 

I  f  ,  during  the  calculation  of  thermal  stress,  we  calculate  the  nosetip 
contour  at  different  points  or  origins,  we  must  use  the  Finite  Element  Method 
to  perform  divisions  -for  all  of  the  prior  mesh  points,  and  this  is  indeed 
quite  cumbersome.  Even  if  we  speak  in  terms  of  a  single  origin  for  the 
nosetip  contour,  because  of  ablation  and  the  various  conditions  that  are 
subject  to  change  at  any  time,  we  still  must  be  able  to  figure  the  origin  of 
the  nosetip  contour  at  any  instant,  and  all  of  these  calculations  are  done 
using  the  Finite  Element  Method's  long-hand  division,  performed  one  by  one, 
by  man,  for  each  mesh  point.  This  makes  the  whole  process  unfeasible. 

Because  we  can  only  perform  calculations  of  thermal  stress  for  a  single, 
fixed  form  or  .for  a  few  forms,  each  calculation  must  have  its  own  division 
done  one  at  a  time.  As  a  result  of  this,  the  development  of  an  automatic 
mesh  dividing  technique  for  the  physical  contour  changes  that  occur  became  a 
necessity.  During  the  compilation  of  the  comprehensive  analysis  algorithm, 
we  discovered  an  automatic  mesh  dividing  technique.  By  using  this  technique 
to  perform  the  divisions  in  the  Finite  Element  Method,  not  only  can  ablation 
changes  be  calculated  for  the  physical  contour,  but  the  contour  changes  along 
the  physical  surface  can  be  matched  with  these  calculations.  Furthermore,  we 
need  only  specify  the  origins  for  the  dimensions  of  the  physical  contour, 
then  for  different  sized  and  shaped  contours  the  automatic  mesh  dividing 
technique  will  perform  the  calculations  and  match  things  up.  As  expressed  in 
Figure  5,  after  a  suitable  selection  of  the  constants  ,  z  and  Q ^  * 

listed  in  the  figure,  we  can  carry  out  the  calculations  for  each  constant  s 
physical  changes  in  seque  :e  by  using  the  Finite  Element  Method. 


tc=Zcf(-^  +  (Zc-Ro)-t9jS-Re] 

Figure  5.  Possible  constants  P Q .  :c<  z^ .  P ^ ,  and  ^ 

for  contour  of  SA I -model  noset  id  at  initial  stage, 
n  /"  a 

Based  on  necessity,  selection  o-f  these  possible  constants  at  different  stages 
can  greatly  simplify  the  process  of  calculation.  In  the  final  quantitative 
analysis  of  the  effects  of  various  types  of  boundary  conditions  produced, 
this  may  provide  significant  contributions  to  engineering  design. 

In  starting  the  comprehensive  analysis  algorithm  for  calculation  during 
initial  re-entry,  based  on  the  external  constants  provided  for  the  thermal 
environment  and  the  projectile  type,  we  can  begin  to  calculate  ablation  and 
the  instantaneous  temperature  fields  and  the  thermal  environment  as  they  are 
coupled  together  in  real-time  CTNssic].  In  fundamental  calculations  of  the 
instantaneous  contour  for  ablation,  and  for  the  instant  temperature  fields 
along  the  contour  for  changing  boundary  conditions,  this  method  of 
calculation  can  provide  for  the  stress  and  strain  fields  at  any  instant 
becuase  of  the  automation  of  the  mesh  dividing  technique  for  the  Finite 
Element  Method.  We  can  use  the  algorithm  to  correctly  calculate  thermal  flow 
with  the  f  o 1 1 ow i n  g  conditions: 

Due  to  the  calculation  of  temperature  and  contour  at  all  times,  the  time 
intervals  were  set  to  a  .2  second  tolerance,  and  resulting!/,  all  constants 
for  calculations  of  the  projectile  are  under  these  limitations. 
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Cone  1  us i on 


This  paper  iterated  the  comprehensive  calculation  algorithm  -for 
ablation,  thermal  environment,  conduction,  and  thermal'  stress.  Furthermore, 
within  the  comprehensive  analysis  algorithm,  an  automatic  mesh  dividing 
technique  was  tested.  Utilization  o-f  the  comprehensive  analysis  algorithm 
tor  carbon-based  noset ips  during  re-entry  ablation,  conduction,  and  thermal 
stress  at  various  atmospheric  levels  can  be  accomplished  with  comprehensive 
numeric  data  analysis.  Moreover,  the  form  and  deterioration  o-f  the  noset  ip 
can  be  calculated  in  succession  using  selective  calculation  techniques  CTN: 
sic]. 
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